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The BRST quartet mechanism in infrared Landau gauge QCD is investigated. Based on the observed positivity violation for 
transverse gluons At, the field content of the non-perturbative BRST quartet generated by An is derived. To identify the gluon's 
BRST-daughter state as well as the Faddeev-Popov-charge conjugated second parent state, a truncated Bethe-Salpeter equation for 
the gluon-ghost bound state in the adjoint colour representation is derived and studied. This equation is found to be compatible with 
the so-called scaling solutions of functional approaches. Repeating the same construction for quarks instead of At, leads in a similar 
way to a truncated Bethe-Salpeter equation for the quark-ghost bound state in the fundamental representation. Within the scaling 
solution the infrared divergence of the quark-gluon vertex is exactly the right one to make this Bethe-Salpeter equation infrared 
consistent. 
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I , As it is well-known by now, gauge-fixing a Yang-Mills (YM) 
Ontheory leaves nevertheless an important invariance: The YM ac- 
^ tion is invariant under a transformation which may be pictured 
i_i as a gauge transformation but with the Faddeev-Popov ghost as 
transformation "parameter" lUl IH [sl- This BRST transforma- 
tion is on the perturbative level the corner stone in the proof of 
renormalizibility, for an accurate discussion see e.g. Ref. ||4j]. 
In addition, it can be employed to define a physical state space 
which in this language is identified with the cohomology of the 
BRST charge operator, see Ref. [|5i] and references therein. Ob- 
PsJ viously, the definition of a positive-definite state space is inti- 
mately related to confinement: The coloured states underlying 
confinement cannot be attributed to asymptotic states in the S - 
. . matrix. To this end it is interesting to note that the BRST coho- 
mology automatically contains only colour-singlet states. 

These remarks make evident that the two purposes BRST 
cohomology serves require two different kinds of BRST mul- 
tipletsQ There are on the one hand the perturbative BRST mul- 
tiplets, first of all the elementary BRST quartet [5] which takes 
care of the cancellation of longitudinal and time-like gluons as 
well as ghosts and antighosts in all physical states. This is noth- 
ing else than the generalization of the Gupta-Bleuler formalism 
of QED. Or, as Weinberg L4J phrases it, the perturbative BRST 
quartet mechanism in YM is an (N^ - l)-fol(fl duplication of 
the single cancellation mechanism in QED. For the definition 
of the state space as BRST cohomology these considerations 
are by far not sufficient. On the contrary, in QCD we expect 
the BRST cohomology to contain only glueballs and hadrons. 



X 



Besides the colourless BRST singlets there are only the so-called BRST 
quartets, see below. 

-N^ - 1 being the dimension of the adjoint representation. 



i.e. bound states. As bound states occur only beyond pertur- 
bation theory it becomes evident that one has to employ non- 
perturbative techniques in studying the BRST cohomology in 
general. 

At this point it is interesting to note that recently a possibil- 
ity has been suggested [6J to avoid the Neuberger 0/0 problem 
of lattice BRST iQ] in a modified lattice Landau gauge which 
employs stereographic projections. Besides other fundamental 
issues elucidated in ref. |6] the following aspect is important 
in the present context: It is now clarified that the global BRST 
charge, necessary for the perturbative as well as for the non- 
perturbative quartet mechanism, can be well defined. 

One interesting aspect of BRST multiplets follows from the 
nilpotency of the BRST transformation, see Sect. 2 for details. 
Every non-singlet state can then produce only one further state, 
making thus a doublet. It proves, however, useful to form quar- 
tets such that the Faddeev-Popov charge conjugated state of the 
daughter state in this doublet is used as a parent state which un- 
der BRST generates the 2nd daughter and thus completes the 
quartet. In this letter we will argue that positivity violation of 
transverse gluons imply that they are parent states. As we will 
show this implies that the 1st daughter has to be a gluon-ghost 
bound state, and the 2nd parent accordingly a gluon-antighost 
bound state. In Landau gauge which is ghost-antighost symmet- 
ric these two bound states become degenerate and are described 
by the same equation. Formally, the same remarks apply for 
the quartet generated by quarks. It is one of the main results 
to be presented here that the structures of the corresponding 
Bethe-Salpeter equations (truncated to keep the infrared lead- 
ing terms) become very similar due to the dressing of the quark- 
gluon vertex. 

In Sect. 2 we demonstrate briefly how to obtain the field 
content of the gluon- and the quark-generated non-perturbative 
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BRST quartetsll In Sect. 3 we discuss the ghost-gluon Bethe- 
Salpeter equation and provide arguments for a massless ghost- 
gluon bound state. In Sect. 4 the ghost-quark bound state equa- 
tion is given and discussed. In Sect. 5 we conclude and close 
with an outlook. 

2. BRST quartets in Landau gauge QCD 

Throughout this letter we will only consider Landau gauge. 
To emphasize the nilpotency of the BRST transformation we 
will work in a representation with Nakanishi-Lautrup field B" 
which becomes on-shell identical to the gauge fixing condition, 
B" - (1/^) c?/jA^ where ^ is the gauge parameter of linear co- 
variant gauges. 

It is useful to picture the BRST transformation as a "gauge 
transformation" with a constant ghost field as parameter 

SbAI = %Dfc^J , 6Bq = -igfZ, c"qA, 

Sbc" = - f/"'"^Zi ci^c' A , 6bc" ^ B"A, (1) 
SbB" = 0, 

where Z)°* is the covariant derivative. The parameter A lives in 
the Grassmann algebra of the ghost fields c", it carries ghost 
number Npp = - 1 . Zi and Z3 are the ghost-gluon-vertex and 
the ghost wave function renormaUzation constants. In Landau 
gauge one has Z] = 1 . 

Via the Noether theorem one may define a BRST charge op- 
erator Qb which in turn generates a ghost number graded alge- 
bra on the fields, 6b^ = {iQe, 't)- Defining the ghost number 
operator Qc one obtains 

Ql^O, [iQ,, Qb] = Qb ■ (2) 

This algebra is complete in the full (indefinite metric) state 
space of YM theory, resp., QCD. The BRST cohomology is 
then constructed as follows: The semi-definite physical sub- 
space KqtQb is defined on the basis of this algebra by those 
states which are annihilated by the BRST charge Qb, Qb\4') - 
0. Since 2| = 0, this subspace contains the space Im Qb of the 
so-called daughter states Qb\(P) which are images of their par- 
ent states in full state space. A physical Hilbert space is then 
obtained as the space of equivalence classes; 

>H(eB) = KereB/ImeB. (3) 

This Hilbert space is isomorphic to the space of BRST singlets. 
All states are either BRST singlets or belong to quartets, this ex- 
hausts all possibilities. Note that the condition QbVIj) - elim- 
inates half of these metric partners from all S -matrix elements, 
leaving unpaired states of zero norm which do not contribute to 
any observable. 



^The main results of this section have been known, of course, since quite 
some time. However, in the way needed in the following it is not easily acces- 
sible in the literature. A more detailed description of the material summarized 
in Sect. 2 may be found in Ref. (31 ■ 



For constructing the so-called elementary quartet [9, 5] one 
considers the asymptotic states related to the time-like and lon- 
gitudinal gluons as well as the ghost and the antighost0 Hereby 
one gluon polarization and the antighost provide the parent 
states, the orthogonal gluon polarization and the ghost yield 
the daughter states. In all physical states the contribution of 
this quartet cancels similar to the cancellation of time-like and 
longitudinal photons in QED. This quartet is strictly perturba- 
tive in the sense that it also exists in the limit g — > 0. From 
the construction of this elementary quartet one can infer how 
to construct other perturbative "multi-particle" BRST quartets 
Starting from a state with negative norm (1st parent) one 
obtains the 1st daughter by acting with the BRST charge opera- 
tor on the 1st parent. The Faddeev-Popov charge reflected state 
of the 1 St daughter provides the 2nd parent. Acting on it with 
the BRST charge operator provides the 2nd daughter 

At the perturbative level the transverse gluons belong to the 
BRST cohomology which is, of course, in open conflict with 
the observed confinement of gluons. Therefore it has been spec- 
ulated already decades ago that the transverse gluons are also 
part of a BRST quartet [9] which is then in turn believed to be 
an important aspect of gluon confinement. At approximately 
the same time it has been observed flUl that the antiscreen- 
ing of gluons (which is a very welcome property as it explains 
asymptotic freedom) is already at the perturbative level in con- 
flict with the positivity of the gluon spectral function. By now 
there is no doubt any more that the transverse gluons of Lan- 
dau gauge QCD are positivity violating, see e.g. Ref. [12] and 
references therein. With the remarks given above this makes 
plain that "one-transverse-gluon" states are BRST parent states. 
Their respective daughters, however, cannot be the elementary 
"one-ghost" states because these are members of the elemen- 
tary quartet. On inspecting eq. ([T]i one can immediately con- 
clude that the corresponding daughter state needs to have the 
field content _ 

Zxf'^'A'^c'' . (4) 

As for every "one-transverse-gluon" state there should occur 
exactly one daughter state this requires the existence of a ghost- 
gluon bound state in the adjoint representation. In this sense the 
resulting BRST quartet is strictly non-perturbative as the forma- 
tion of bound states can be described only with non-perturbative 
techniques. The Faddeev-Popov charge reflected state is then 
an antighost-gluon bound state. Here Landau gauge provides 
an advantage as compared to general linear covariant gauges: 
In the limit ^ ^ the formalism becomes ghost-antighost- 
symmetric, and thus the existence of a ghost-gluon bound state 
implies the occurrence of a degenerate antighost-gluon bound 
state with same quantum numbers. Even having then the 2nd 
parent, the BRST transformation ([TJ leaves then three possi- 
bilities for the 2nd daughter: A ghost-antighost bound state, a 
ghost-antighost-gluon bound state, or a bound state of two dif- 
ferently polarized gluons. However, studying the 2nd daughter 



As the Nakanishi-Lautrup field 6" relates to a linear combination of time- 
like and longitudinal gluons, the so-called backward polarization, the corre- 
sponding asymptotic "one-gluon" parent (daughter) states are forwardly (back- 
wardly) polarized, see e.g. Chapter 16 of Ref. 1 10]. 
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state is beyond the scope of this letter, and we will return below 
to the ghost-gluon bound state. 

The respective issues for quarks in the Landau gauge are 
much less clear. First of all, it is not known whether quarks 
violate positivity. Although for light quarks dynamical chi- 
ral symmetry breaking (and for heavy quarks explicit chiral 
symmetry breaking) determines the infrared behaviour of the 
quark propagator the analytic structure of the quark propaga- 
tor is highly sensitive to details in the quark-gluon vertex, see, 
e.g., Ref. il3ll . The quark-gluon vertex for light quarks is, on 
the other hand, also very strongly influenced by dynamical chi- 
ral symmetry breaking 15, 16]. Second, an inspection of 
eq. O reveals that if a BRST quartet is generated by quarks 
it can only be a non-perturbative one. Which mechanism then 
guarantees that the corresponding bound states are degenerate 
with the quark states is completely unknown. 

Before trying to give an at least partial answer to this we will 
discuss the infrared behaviour of Landau gauge YM theory and 
the implications for a possible ghost-gluon bound state and its 
role in the corresponding BRST quartet. 

3. Landau gauge YM theory and ghost-gluon bound state 

Over the last decade the infrared behaviour of Landau gauge 
YM theory has been in the focus of many studies. Hereby 
it is interesting to note that in the deep infrared quite general 
statements can be deduced by employing functional equations. 
On the one hand, Dyson-Schwinger equations have been used 
to extend a previous analysis of gluon and ghost pro pag ators 
1 13, HE Sim to aU Yang-Mills vertex functions iHIBIl . 
Employing in addition the Function Renormalization Group 
Equations, and requiring that these two, seemingly diff'erent, 
towers of equations have to provide identical Green's functions, 
allows a powerful restriction on the type of the solution: There 
is one unique scaUng solution with power laws for the Green's 
functions fl^, ^23] and a one-parameter family of solutions, the 
so-called decoupling solutions. These infrared trivial solutions 
possess as an endpoint exactly the scaling solution which is 
characterized by infrared power laws. Numerical solutions of 
the decoupling type (there called "massive solution") have been 
published in 11261 12711 and references therein. A recent detailed 
description and comparison of these two types of solutions has 
been given in Ref. [28 1, see also Refs. Ii29..30, 3L.321 . 

Although almost all lattice calculations of the gluon propa- 
gator favor a decoupling solution lattice studies at strong cou- 
pling 11331 134 13511 reveal the existence of a regime where the 
scaling relation between the gluon and the ghost propagators 
is fulfilled, and the corresponding infrared exponent k is very 
close to the value determined in truncated continuum studies 
with K = 0.595. A potential explanation of these conflicting re- 
sults has been discussed recently ^36*]: The infrared behaviour 
of the Green's functions may depend on the non-perturbative 
completion of the gauge. 

In the present context it is important to note that the scaling 
solution respects BRST symmetry whereas every decoupling 
solution breaks it ll28[l . although very likely only softly. In a 
strict sense our analysis as presented in the following will be 



only valid if the scaling solution is a correct one. However, the 
situation is not as drastic as it may seem. First, if the conjec- 
ture of Ref. [36] is correct it is sufficient that only one non- 
perturbative completion of Landau gauge with scaling solution 
exists to make our analysis well-founded. Second, even if only 
decoupling type of solutions will turn out to be correct an ex- 
tended BRST-like nilpotent symmetry is likely to take the role 
of the BRST symmetry 114 ill , or the soft BRST symmetry break- 
ing can be treated as spontaneous symmetry breaking [42], see 
also Ref. Ii43i] and references therein. Note also that all our ar- 
guments about infrared dominance of diagrams will likely stay 
correct: The numerical value of a diagram which is infrared 
leading in the scaling solution will be enhanced {i.e. numeri- 
cally large) in a decoupling solution if the latter is not far from 
the endpoint given by the scaling solution. 

As a basis for the following discussion we will summarise 
the infrared behaviour of all one-particle irreducible Green's 
functions in the scaling solution in the simplified case with only 
one external spacelike scale — > 0. For a function with n 
external ghost and antighost as well as m gluon legs one has: 



2\(n-m)K 



(5) 



This solution fulfills all functional equations and all Slavnov- 
Taylor identities. It verifies the hypothesis of infrared ghost 
dominance (yf] and leads to an infrared diverging ghost propa- 
gator as well as infrared diverging three- and four-gluon vertex 
functions. 



As already stated gluons violate positivity 11121 11311 . For the 
scaling solution this can be immediately deduced from the fact 
that for this solution the gluon propagator vanishes at zero vir- 
tuality, p- - 0. A further important property of the scaling so- 
lution is the infrared trivial behaviour of the ghost-gluon vertex 
which is in agreement with general arguments j 2^ . 

To search for the anticipated ghost-gluon bound state we 
want to truncate the ghost-ghost-gluon-gluon scattering ker- 
nel to the infrared leading term. To this end we employ first 
the MATHEMATICA package DoDSE, resp., DoFun 139", '4^ 
to derive the diagrammatic expressions for the one-particle 
irreducible ghost-ghost-gluon-gluon scattering kernel. The 
diagram-by-diagram infrared power counting^ verifies that in 
the scaling solution the infrared exponent is -k. It also provides 
the infrared leading terms. 

As two different kind of fields are involved there exists two 
distinct possibilities for the Dyson-Schwinger equation accord- 
ing to which leg one puts the bare vertex. Placing the bare ver- 
tex on a gluon line gives 56 diagrams (59 if quarks are included) 
with six (seven) of them containing 5-point or 6-point functions. 
Eleven of them are infrared leading. Assigning the bare vertex 
to a ghost leg leads to a Dyson-Schwinger equation with 1 3 dia- 
grams on the r.h.s, one containing a 5-point function, and all of 
them infrared leading. As we have to truncate the system such 
to neglect n > 5-point functions it is obvious that the second 
choice minimizes the truncation errorj^ This is substantiated by 



""Some more details of this analysis are given in ref. |8|. 

^Note that either neglecting or approximating the one-particle irreducible 5- 
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Figure 1: Graphical representation of tlie gluon-ghost Bethe-Salpeter equation. A cross denotes a dressed vertex. 



the fact that the first choice would provide eventually a Bethe- 
Salpeter equation with one less term as if we will employ the 
second choice, see below. 

The requirements for the diagrams on the rh.s to be kept 
are: It should contain the one-particle irreducible ghost-ghost- 
gluon-gluon four-point function and no « > 5-point function, it 
should be infrared leading, and it should contain the interaction 
in the ghost-gluon channel. This leaves two diagrams: (i) one 
with two ghost and one gluon propagator on internal lines. This 
is effectively a ghost exchange. (NB: A more precise descrip- 
tion is that the leading tree-level diagram describes the splitting 
of the incoming gluon into a ghost-antighost pair and a fusion of 
the incoming ghost with the exchanged (anti-)ghost to a gluon, 
cf. the upper rh.s of Fig. [T]) (ii) one with two gluon and one 
ghost propagators on internal lines. This is a gluon exchange. 
Note that this diagram is infrared leading because in the scaling 
solution the fully dressed three-gluon vertex is infrared diver- 
gent with an exponent -3/c If we had chosen to put the bare 
vertex on the gluon leg (cf. the first choice above) this diagram 
would be infrared suppressed. We would also like to remark 
that from the 14 tensor structures of the three-gluon vertex at 



point function is necessary to obtain a solvable Bethe-Salpeter equation. How- 
ever, as the diagram containing the 5-point function possesses the same infrared 
exponent as the other diagrams the resulting equation stays qualitatively correct, 
and it is reasonable to assume that it is quantitatively satisfactoiy. Please note 
also that the truncation is not based on the number of diagrams but based on the 
argument to keep as many infrared leading diagrams as technically possible. 



least ten are infrared divergent with the same exponent ||44 1 . 

Assuming the existence of a bound state as well as employing 
the usual decomposition of the (ghost-ghost-gluon-gluon) four- 
point function into Bethe-Salpeter amplitudes and performing 
the expansion around the pole (see e.g. Sect. 6.1 of Ref. |45]) 
one arrives at the Bethe-Salpeter equation depicted in Fig. [1] 
Using the propagator parameterizations of e.g. Ref. |13], the 
ghost-gluon vertex of Ref. ll46tl . and the three-gluon vertex of 
Ref. |44] one can derive a self-consistent equation for the corre- 
sponding Bethe-Salpeter amplitude containing otherwise only 
known quantities. This leads, mostly due to the many tensor 
terms of the three-gluon vertex, to very lengthy expressions 
which will be given elsewhere. For the illustrational purposes 
of this letter we will neglect the gluon-exchange term and keep 
only the ghost exchange. 

As stated above in Landau gauge the ghost-gluon vertex stays 
even in the scaling solution infrared trivial. Thus it is suffi- 
cient to consider only bare ghost-gluon vertices, i.e. it is suf- 
ficient to consider the ladder approximation to the first term of 
the above Bethe-Salpeter equation. (We also use Zi = 1 llssll .) 
The explicit expression for the kernel, assuming as gauge group 
SU(A^c) and working in Euclidean momentum space, reads then 

Daiqi - k2)igff'q2f, (6) 
= -g^Nl5''''r,q2,Dc{r\ 
where we have already taken care of the re-projection onto the 
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adjoint representation and have used the momentum assignment 
as in the rh.s of Fig. [1] Hereby Ddqi -^2) - Dc(r) is the ghost 
propagator 

With this kernel and denoting the gluon propagator as D^y{k) 
we arrive at the gluon-ghost Bethe-Salpeter equation for bound 
state with four-momentum P 

—I-r,q2^Dc{r) (7) 
D,Ak2)Dc(q2W^p(k2,q2;P). 

As a side remark we want to mention that to arrive at the equa- 
tion for the antighost-gluon bound state we only need to invert 
all ghost momenta and obtain therefore: 



Tf,p(ki,qi;P) 



(r) 



(8) 



Dy^(k2)Dc(q2Wa-p(k2, q2, P)- 

As argued above we are looking for a massless bound state. 
This allows to specialise to the soft limit P ^ 0. Employing 
the transversality of the Landau gauge gluon propagator, 

/ kykAzie) 

DyAk) = \6y^ - -j^ j -j^ (9) 

it is straightforward to show that T^^pik, -k; 0) (resp., r^p(A:, k; 0) 
is transverse, too: 



r,p(k,-k;0) 



'MP 



'-^\F(k\ 



(10) 



Defining the ghost renormalization function G(q^) - q^Dc(q^) 
one obtains then from eq. dTji 



^(^1 ■ k2) 



k2 G((ki + k2f) G(kl) Z(kl) 



{ki- 
(ki 



■k2)^ 

■k2f^ 



k^ 
''2 

Fikf). 



k^ 
'^2 



(11) 



It is an easy exercise to verify that the Bethe-Salpeter equation 
([8]l provides exactly the same equation as ( fTTT i and therefore the 
degeneracy of the massless antighost-gluon bound state if the 
massless gluon-ghost bound state exists. 

The numerical solution of eq. (fTTl i is beyond the scope of 
the present letter, it will be presented elsewhere. Here we only 
want to mention that, first, F(k^) has the peculiar feature F(Q) = 
0, and, second, due to truncation errors we do not expect the 
equation to be exactly fulfilled. As, however, we have kept the 
infrared leading term it is important to test whether there is no 
contradiction on the analytical level. 

For the scaling solution we have the power laws 



Zip') 



(p'r. 



(12) 



which implies that the kernel of eq. (fTTT i possesses a vanishing 
anomalous infrared exponent and is independent of the infrared 
exponent k. We therefore get the expected result that the Bethe- 
Salpeter equation for a massless (anti-)ghost-gluon bound state 



is consistent within the infrared analysis of the scaling solu- 
tion of functional equations. Furthermore, as the combination 
a^^(p^) = j^G^{p^)Z{p^) has the form of a running coupling 
with an infrared fixed point a^''{Q) - 8.92/A^c (see e.g. Sect. 2.3 
of Ref. II47II ) it is plausible that the kernel of eq. (fTTT i is strong 
enough to lead to a bound state. All these findings corroborates 
the validity of the employed method, and therefore we will in- 
vestigate the quark-gluon bound state with the same method. 

4. Quark-gluon bound state 

The scaling solution for the YM Green's functions leads to 
dynamical chiral symmetry breaking in the quark sector LI61I . 
The quark propagator which is of the general form 



S (P) = ..7 . ..O...O. ^f(P ) 



(13) 



p2 + mHp^) 

is then infrared finite, 5(0) - Z/(0)/M(0). The twelve possible 
Dirac tensor structures of the quark-gluon vertex are then all 
infrared divergent with an infrared exponent -k - 1/2. As a 
side remark we want to add that the same infrared divergence 
results for vanishing gluon momentum, and that this leads to 
l/k'^ behaviour of the kernel in the four-quark function, k being 
the momentum exchange. This is indicative of a linearly rising 
potential between static quarks, and thus quark confinement. 

Furthermore, the Slavnov-Taylor identities require that the 
ghost-ghost-quark-quark scattering kernel is infrared trivial, see 
Sect. 3.9 in Ref. [16]. This information is absolutely necessary 
for the following analysis. 

As in the previous subsection we have two choices for the 
Dyson-Schwinger equation according to which leg one puts the 
bare vertex. Both choices lead to seven diagrams on the rh.s 
of the equation, in both choices there appears one diagram with 
a 5-point function (ghost-ghost-quark-quark-gluon). However, 
choosing the quark leg to carry the bare vertex only one dia- 
gram is infrared leading, namely exactly the one with the 5- 
point function. Thus truncating the equation to n < 4-point 
functions would lead to a contradiction. Choosing a ghost leg 
to place the bare vertex there are four infrared leading diagrams: 
One with the 5-point function, and three according to the three 
possible (s,t,u) interaction channels. Therefore, employing the 
same truncation requirements and the same derivation of the 
Bethe-Salpeter equation as in the previous subsection one ar- 
rives at the equation depicted in Fig.|2] This equation is in full 
agreement with the infrared analysis of the scaling solution, i.e. 
it is a valid bound state equation, and in its kernel the infrared 
exponent k cancels. 

Employing all twelve tensor structures of the quark-gluon 
vertex does generate a very lengthy equation. Again we will 
restrict for illustrational purposes to the perturbatively leading 
tensor component. 



r"^(puP2;k) -ig—y^V(pup2;k) 



(14) 



where V diverges like {k^) " '^^ for k^ — > 0. Using the result 

A" A' A'' 



2 2 2 



j^abc 



1) 



(15) 
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Figure 2: Graphical representation of tlie quark-gliost Betlie-Salpeter equation. Note that the quark-gluon vertex is fully dressed. 



valid for SU(A^c) one can project the bound state onto the funda- 
mental colour representation. The Bethe-Salpeter equation for 
the ghost-quark bound state amplitude with total momentum P 
reads then: 



(16) 



V(pi 



(In) 

P2\ k)S{p2)Dc{qiW{p2, qi; P) 



with p2 = pi + k and q2 - q\ - k. At first sight it may look 
like that the Bethe-Salpeter amplitude of this scalar-fermion 
bound state has four independent amplitudes. This is, however, 
not correct. Should the bound state occur at a non-vanishing 
P^ - Mg < (to avoid saying the "ghost-quark bound state 
is massive") there will exist two degenerate solutions with re- 
spectively positive and negative energy. Then it is sufficient to 
study only the positive-energy solution by explicitly putting a 
positive-energy projector into the tensor decompositioiu 



np,q;P)^{H[ + —{p-4)H', 



A+ 



= ^(1 + 



-). 



(17) 

The Dirac algebra to generate the set of coupled equations for 
//J and //^ is lengthy but straightforward. However, as the equa- 
tions themselves are somewhat vast we refrain from displaying 
them here. 

Again we want to simplify for illustrational purposes and 
simply set Y(p,q;P) - H(p,q,P), i.e. to a scalar function. 
In the case of the dominance of the "upper Dirac component" 
this should be a reasonable approximation. To further demon- 
strate the properties of the kernel we make the simplest pos- 
sible choice for V(pi, p2',k) agreeing with its infrared limit: 
y(Pu P2\k) — > G{k^)l Vfe^. The corresponding equation for the 
simplified amplitude H{p, q; P) reads then: 

,A^3 - 1 r d% Z(k^) G(k^) G{q\) 



Hipuqi;P) = 



1 r d^k 

~ J (2? 



Z/(P2) 

pI + mHpD 



)4 k^ 



27 Kl -Pl 



qi ■ p\p\ ■ P2^ 



k^ 



H(p2,q2;P). 



(18) 



'This is very similar to the bound state equation of a scalar diquark and a 
quark for a nucleon, see e.g. Sect. 7.3 of Ref. |45|| . 



The Structural similarities to eq. (ITTl are striking. Again it is 
likely that the kernel is strong enough to lead to a bound state. 
Providing further evidence that these bound states are degen- 
erate with quark states will, however, require the study of the 
system with the full ansatz ( fTTl l. 

At the end of this section we want to add a cautionary remark: 
Although the statements given in this section do not hold only 
in the quenched limit of QCD but also for QCD with dynamical 
fermions, the resulting equations cannot be used directly. Note 
that e.g. in the step from eq. (fTSI l to eq. (fTSl l a result for the 
quark-gluon vertex lfl6ll has been used which is only known in 
the quenched limit. 

5. Conclusions and outlook 

The observed positivity violation of the Landau gauge gluon 
propagator imply that, if BRST is an unbroken symmetry, the 
transverse gluons generate a non-perturbative BRST quartet. 
The other three members are bound states. For two of these 
we have derived a (truncated) bound state equation by analysing 
the one-particle-irreducible ghost-ghost-gluon-gluon four-point 
function. The corresponding Bethe-Salpeter equation contain- 
ing in its kernel infrared divergent Green's function is infrared 
consistent, a result which is far from being trivial. In addition, 
we have provided arguments that the kernel is strong enough 
such that a bound state is formed. All these findings corrob- 
orate the existence of the non-perturbative BRST quartet with 
transverse gluons as first parent states. On the other hand, one 
might also argue that the employed method is consistent and 
can therefore be extended to study the analogue problem for 
quarks. 

With the same method we derive then the bound state equa- 
tion for a quark-ghost bound state truncated to the infrared lead- 
ing term. We note that this equation is within the scaling solu- 
tion only then infrared consistent if the recently found infrared 
divergence of the fully dressed quark-gluon vertex is taken into 
account. We take this as further evidence that the infrared be- 
haviour of the quark-gluon vertex is intimately related to the 
issue of quark confinement. Nevertheless, further studies of 
the quark-gluon vertex are highly desirable to generate progress 
on the question whether in Landau gauge QCD there is a non- 
perturbative BRST quartet generated by quarks. 
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The numerical solution of the gluon-ghost bound state equa- 
tion (fTTT l can be done with standard techniques, it will be pub- 
lished elsewhere. In order to solve numerically the quark- 
ghost bound state equation ( fTSI l and its analogue for the quark- 
antighost bound state the full Dirac tensor ansatz (fTTl l has to 
be plugged in, and the corresponding Dirac algebra has to be 
performed. However, although the resulting set of equations 
are lengthy they can be solved with methods which are by now 
standard. 

Here we did not touch on the question of the respective sec- 
ond BRST daughters. Following the same arguments as em- 
ployed here the corresponding tasks should be doable. How- 
ever, before starting an investigation in this direction it might 
be worthwhile to repeat the presented investigation with Func- 
tional Renormalization Group Equations to minimize further or 
even eliminate truncation errors. 
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